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Abstract 



We examine the contributions of soft gluons to the Higgs production cross 
section at the LHC in the Standard Model and its minimal supersymmetric 
extension. The soft gluon radiation effects of this reaction share many features 
with the Drell-Yan process, but arise at lowest order from a purely gluonic initial 
state. We provide an extension of the conventional soft gluon resummation 
formalism to include a new class of contributions which we argue to be universal, 
and resum these and the usual Sudakov effects to all orders. The effect of these 
new terms is striking: only if they are included, does the expansion of the 
resummed cross section to next-to-leading order reproduce the exact result 
to within a few percent for the full range of Higgs boson masses. We use 
our resummed cross section to derive next-to-next-to-leading order results, and 
their scale dependence. Moreover, we demonstrate the importance of including 
the novel contributions in the resummed Drell-Yan process. 
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1. Introduction 



The search for Higgs particles |l| is one of the most important endeavors for 
future high energy e + e~ and hadron collider experiments. The Higgs boson is the 
only particle of the Standard Model (SM) which has not been discovered so far. The 
direct search at the LEP1 experiments via the process e + e~ — > Z*H yields a lower 
bound on the Higgs mass of 65.2 GeV f2|. Theoretical consistency restricts the Higgs 
mass to be smaller than ~ 700 GeV 0. The dominant Higgs production mechanism 
at the LHC, a pp collider with a cm. energy of 14 TeV, is the gluon fusion process 
gg — > H which is mediated by a heavy quark triangle loop at lowest order [|J. As 
an important step to increase the theoretical precision the two-loop QCD corrections 
have been calculated, resulting in a significant increase of the predicted total cross 
section by about 50 - 100% f5|||. The dependence on the unphysical renormalization 
and factorization scales decreased considerably by including these next-to-leading- 
order (NLO) corrections, resulting in an estimate of about 15% for the remaining 
scale sensitivity ||. It is important to note, and we will demonstrate, that the NLO 
corrections are dominated by soft gluon radiation effects. 

The minimal supersymmetric extension of the Standard Model (MSSM) is among 
the most attractive extensions of the SM. It requires the introduction of two Higgs 
doublets leading to the existence of five scalar Higgs particles, two scalar CP-even 
h, H, one pseudoscalar CP-odd A and two charged bosons H ± . This Higgs sector can 
be described by fixing two parameters, which are usually chosen to be tg/9, the ratio of 
the two vacuum expectation values, and the pseudoscalar Higgs mass Ma- Including 
higher order corrections to the Higgs masses and couplings up to the two-loop level, 
the mass of the lightest scalar Higgs particle h is restricted to be smaller than ~ 130 
GeV [0. The direct search at LEP1 sets lower bounds of about 45 GeV on the masses 
of the MSSM Higgs bosons |2|. The dominant neutral Higgs production mechanisms 
at the LHC are the gluon fusion processes gg — > h,H,A and the associated production 
with a bb pair gg, qq — > bbh, bbH, bbA which becomes important only for large tg/3 [|J. 
The coupling of the neutral Higgs particles to gluons is again mediated by top and 
bottom loops, with the latter providing the dominant contribution for large tg/5, and 
squark loops, if their masses are smaller than about 400 GeV |J. (In this paper 
we shall assume the squark masses to be 1 TeV, so that squark loops can safely be 
neglected.) The two-loop (NLO) QCD corrections to the gluon fusion mechanism 
have also been calculated [§,|9| and conclusions completely analogous to the SM case 
emerge. Soft gluon radiation effects again provide the dominant contribution to these 
corrections, for small tg/3. For large tg/3 bottom mass effects will be of similar size 
due to the dominance of the bottom quark loops and we will not consider this regime 
in this paper. 

In previous analyses the resummation of soft gluon radiation in the transverse 
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momentum distribution of the Higgs bosons has been performed [|Kj, which are sig- 
nificant at small p?. We consider universal soft gluon effects on the total production 
cross section and demonstrate that these dominate the NLO corrections both in the 
SM and the MSSM for small tg/3. The study of these effects in higher orders, and 
their resummation to all orders is our purpose in this paper. 

Although our main focus is Higgs production, we will consider soft gluon effects 
in the Drell-Yan process for comparison. Higgs production shares many features 
with this reaction, apart from the species of leading initial state partons, e.g. it 
also proceeds at lowest order via a color singlet hard scattering process, and is a 
2 — > 1 process at lowest order. The Drell-Yan process has been studied by performing 
exact perturbative QCD calculations up to next-to-next-to-leading order (NNLO) in 



Ref. [11] and in the context of soft gluon resummation in [12|. In this paper we 



present an extension of the conventional soft gluon resummation formalism, in which 
we use the Drell-Yan reaction to gauge its quality and importance. We then apply the 
extension to Higgs production to derive the first estimates of NNLO effects. These 
estimates are important in view of the size of the NLO corrections. 

The paper is organized as follows: In section 2 we describe the construction of the 
resummed cross section and the extension of the soft gluon resummation formalism. 
In section 3 we present NLO and NNLO results from the expansion of the resummed 
cross sections for Higgs production and the Drell-Yan process. We conclude and 
present an outlook in section 4. 



2. The Resummed Exponent 

In this section we derive the resummed partonic cross section for Higgs boson 
production via gluon fusion. In order to set the stage, we must first discuss some 
preliminary approximations to the exact NLO calculation in Ref. 

In the Standard Model, the leading order (LO) process consists of gluon fusion into 
a Higgs boson via a heavy quark triangle loop, see Fig. [l|. Because the Higgs coupling 
to fermions is proportional to the fermion mass, the top quark strongly dominates 
this coupling, constituting about 90% of the total coupling. Our first approximation 
is to neglect the quark-antiquark and quark-gluon channels as these contribute at 
next-to-leading order (NLO) to the full cross section with less than 10% []5J. 

The exact NLO calculation in || was performed for the general massive case, 
i.e. all masses were taken into account explicitly. However, the following very useful 
approximation was identified, involving the heavy top mass limit of the calculation. 
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Figure 1: Higgs boson production via gluon fusion mediated by top- and bottom quark 
loops 



Let us define the (NLO) K-factor^ by 

jst+b / ^ x _ 0~NLo{T t ,T b ) 

K NLoi r t,r b ) = — (1) 

°Lo\Tt, n) 

where CLo/NLoijt,^) denotes the hadronic gluon-fusion cross section for the general 
massive case, calculated exactly in LO/NLO, and the scaling variables are defined by 
T Q = Am%/Ml (Q = t,b). 

In Figs. and |3| we compare aNLoijt, T b) with the approximation 

k nlo(°°) X °~Lo(T t ,T b ) (2) 

for scalar and pseudoscalar Higgs boson production, where the K-factor i^ LO (oo) 
takes into account the top quark contribution to the relative QCD corrections only, 
in the limit of a heavy top quark. We observe that the approximation ([J) is accurate 
to within 10% for the full Higgs mass range Mu > 65 GeV of the SM Higgs boson 
as well as the pseudoscalar Higgs particle of the MSSM for small tg/3 [At the ti 
threshold Ma = 2m t the pseudoscalar cross section develops a Coulomb singularity 
so that perturbation theory is not valid in a small margin around this value for the 
pseudoscalar mass Ma || .] . The same accuracy of the approximation emerges for the 
two scalar Higgs particles of the MSSM for small tg/3. Our second approximation then 
consists of assuming that K f (oo) x atoijt-, T b) is a valid approximation to a{jt, r&) for 
all orders, i.e. we will assume that the higher-order K-factor, when computed in the 
infinite top mass limit and combined with the massive LO cross section, will give a 
good approximation to the higher order cross section in the general massive case. In 
fact, we will see that at NLO the bulk of the K-factor is due to soft and collinear 
gluons, which do not resolve the effective coupling. The assumption that this persists 



2 It should be noted that we use here NLO parton densities and strong coupling a s in the NLO 
cross sections and LO quantities in the LO cross sections for this "hadronic" K-factor. This leads 
to K-factors smaller than two in contrast with the "partonic" K-factor, which is defined with NLO 
parton densities and strong coupling also in the LO cross section. 
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Figure 2: Exact and approximate results in the heavy top quark limit for the total SM 
Higgs production cross sections as a function of the Higgs mass Mh- The top mass 
has been chosen as m t = 175 GeV and the bottom mass as nib = 5 GeV. CTEQ4M 
parton densities jTBtf with NLO strong coupling [a s {M'^) = 0.116/ have been used. 
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Figure 3: As in Fig.^ but now for the MSSM pseudoscalar Higgs particle for three 
values of tg/3. 
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to higher orders is supported by the validity of the infinite top mass approximation 
at NLO. 

In the MSSM, the validity of these approximations depends strongly on the pa- 
rameter tg/3. For tg/3 < 1.6 the top quark contribution to the cross sections amounts 
to more than 70%. The heavy top quark limit is thus a reasonable approximation in 
this regime. For large values of tg/3 the bottom loop contribution becomes significant 
so that the approximations are no longer valid. Still, the infinite top mass approx- 
imation deviates from the full NLO result, including bottom contributions, by less 
than 25% for tg/3 < 5 as can be inferred from Fig.[| 

We are now ready for the construction of the resummed partonic cross section, 
for which we will employ the methods of Ref. ||14|| . In order to retain similarity to the 



Drell-Yan case, we will denote the Higgs mass squared with Q 2 throughout the text 
of this paper. 

In the approximations outlined above, the regularized total Higgs production cross 
sections may be written in d = 4 — 2e dimensions as [(f) = h, H, A] 

^(t&Q 2 , /J 2 ) = f dxi ( dx 2 g{xi) g(x 2 ) &g g {z, Q 2 ,At 2 ,e) (3) 

or, in terms of moments, 

a*(N,Q 2 ,^ 2 ) = £dr* t^ 1 a*{r ch Q\^)=g\N + l) &$ g {N, Q 2 ,fx 2 ,e) (4) 

where = Q 2 /S, z = t^/(xiX 2 ), S is the hadronic cm. energy squared, fi is the 
dimensional regularization scale, and g(x) is the bare gluon distribution function. 
Note that the definition is such that the dependence on the moment variable in the 
parton densities is shifted by one unit compared to the dependence of the partonic 
cross section. In this way we remove an overall 1/z factor, emphasizing the soft gluon 
contribution to the partonic cross section. 

In the approximations discussed in the beginning of this section, the c?-dimensional 
partonic cross sections can be cast into the form 



°tq = a o H PA Z ^ Q 2 /^ 2 , e) (5) 



with the coefficients 



h,H h,H G Fa 2 sB NcC F T 2 ^ ., ,.. 
a ° = 9t 1152^ 7— I") , (6) 
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G F al B N c C F T 2 (l , ,, ... 
= 9t — ~ 1 — 2 ' (') 



where a St B is the bare strong coupling constant (with dimension 2e) and gf(<f) = 
h, H, A) denote the modified top Yukawa couplings normalized to the SM coupling, 
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which are given in ||. The factor in eq. @ stems from the effective coupling of 
the Higgs boson to gluons in the heavy top quark limit, which can be obtained by 
means of low energy theorems |||l5j . They are givenQ by |5l,|l6| 



3tt 



Kh,H 



K A 



a 



(5) (m 2 ) - 1+7»K ; (WI?)] Ws°'(Ml H )J fiatf'W)] 



(5), 



(9) 



where ai 1 ^ (n/ = 5,6) is the strong coupling constant in the MS scheme including 
nf flavors in the evolution; the couplings for 5 and 6 flavors are related by |[Tj| 

(10) 



af\m$ 




/3(a s ) denotes the QCD (3 function and /3 t (a s ) its top quark contribution^, which is 
given by |18],|20| [n/ = 5 is the number of light flavors] 



371 a ( \ 
-o Pt(ot s ) 



a: 



, 19 a s 6793 - 281 n f a 2 _ „. 



7 m (at g ) is the anomalous mass dimension including 6 flavors, which can be expressed 
as 



7m(« s ) 



2— + 

7T 



ra f + 1 1-4 + Of 

12 18 L } 7 7T 2 v 



(12) 



Using these expansions the effective scalar couplings Kh ; H are given byf] 



1 + 



11 ai 5 )(r// 



, , 3866 - 201 n/ fa 



2 7T 144 \ 

153- 19 n, «f(M^)-a( 5 )K) 




33 - 2n 



f 



7T 



+ 0(a2 



(13) 



where m t denotes the pole mass of the top quark. The scale dependence of the strong 
coupling in the lowest order cross section eq. (|7|) and the factor eq. (^) includes 5 
light flavors, i.e. the top quark is decoupled. In the rest of the paper we identify 



a s ( M 2 ) = «( 5 V) 



(14) 



3 The last two factors in the large bracket originate from the anomalous dimension of the gluon 
operators |p^| . The top mass m* denotes the scale invariant MS mass frit = TOt(mt). We would 
like to thank the authors of Ref. |l7j] for pointing out two errors in our treatment of the strong 
coupling in eq.(^) and the anomalous mass dimension of cq.(|l2|) in an earlier version of this paper. 
The numerical size of these errors is about 0.03% and thus negligible. 

4 The factors k^.h include the top quark contribution at vanishing momentum transfer, which 
differs from the top quark contribution to the MS (3 function by a finite amount at 0{a A s ) JnJ. 

5 This result agrees with Ref. p7). 
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Note further that the expression ([5]) is not yet finite for d — > 4; mass factorization 
and renormalization of the bare coupling in the Born cross section will be carried 
out after resummation. In eq. @ we denote the correction factors by p^(z, Q 2 /p 2 , e), 
which are defined in the infinite mass limit without the factorizing corrections to 
the effective coupling. They may be expanded as 



z,Q 2 /p 2 ,a(p 2 ),e) = £ a n (p 2 )pf(z, Q 2 /p 2 , e) 

n=0 



where we define, for the sake of convenience, 

a s (p 2 ) 



a(p ) = 



7T 



(15) 



(16) 



Here ct s (p 2 ) is the renormalized strong coupling and we have chosen the renormal- 
ization scale equal to p for the moment. From Ref. |5||| we can derive the following 
expression for the first two coefficients of the SM Higgs correction factor 

$\z,Q 2 /p 2 ,e) 



p { Hkz,Q 2 /p 2 ,e) 



5(1 -z) 



;i7) 



1 + z 4 + (1 - z) 4 



U+2e 



+ 5(1 



11 203 vr 2 
6e + + y 
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-zHl 



3 " 2e } 



Pa 



C A 5(l-z) 



(19) 



Note that we have implicitly redefined the scale p by p 2 — > p 2 exp[— (ln(47r) —Je)) to 
eliminate factors (47r) e and T(l — e)/T(l — 2e). The plus distribution in eq. (p!8|) is as 
usual defined by 



dzg(z)[f(z)\. 



dz(g(z)-g(l))f(z)-g(l) / dzf(z 



(20) 



We will now construct a resummed expression for p^(z,Q 2 j ' p 2 ,a,e) by means of 
the methods described in Ref. . Near the elastic edge of phase space the Higgs cross 
section in the infinite mass limit may be factorized into hard, soft and jet functions, 
in completely analogy with the Drell-Yan cross section. Following the arguments of 
Ref. [H| this leads to the Sudakov evolution equation 
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P<j>(z, Q 2 /p 2 , a(p 2 ), e)= —W^z', Q 2 /p 2 , a(p 2 ), e)p 4> (z/z', Q 2 /p 2 , a(p 2 ), e) . 

J z Z 

(21) 

In order to solve eq. (21) we must impose a boundary condition. We will shortly 
argue |T3J that we may use the boundary condition p r f ) (N,Q 2 /p 2 = 0,a(p 2 ),e) = 1 
for the moments, or in z-space 



dQ- 



Zy Q 2 /^ 2 = 0, a(p 2 ), e) = 5(1 - z). 



(22) 
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The solution to eq. ( plf) is then, in Mellin space, 



(N,Q*/ t i i ,a(jM i ),e) = exp 



' rQ 2 d£ 2 ~ f t 2 



where as above d = 4 — 2e > 4. The formal Mellin inversion reads 

/too HAT 
-ioo iKl 



(23) 



(24) 



The solution fl23|) may be expressed as 



p^(iV,g7// 2 ,a(/i 2 ),e) 



exp 



fizz 



JV-l 



^Q 2 /m 2 dA 



T {a(A,a(/i 2 ),e)< ) (z,l/i,,e) 



+ « 2 (A, «(/i 2 ), e)Wf l/v, e) + G(a 3 )} 



(25) 



with v{z) an arbitrary function. Note that u(z) is to be treated as part of any 
plus distributions in W^, because it arises from the originally ^-independent scale 
ratio. In deriving eq. ( p5|) we have used the renormalization group invariance of 
the evolution kernel Wq. Note further that we have expanded the full evolution 
kernel Wf(z, £ 2 //U 2 , ct(/i 2 ), e) in the d- dimensional running coupling constant, with 
A = u(z)^ 2 /^ 2 . The defining equation of the d- dimensional running coupling is 

Al _ e 0[A'a(A^a(/z 2 ),6)] = ^ g) _ ^ g) ? (2g) 

with the boundary condition a(l, a(// 2 ), e) = a(/i 2 ). Here 6 2 = (HCU — 2n/)/12 and 
6 3 = 34C1/48 - (20C A /3 + 4C F )n f /32. The solution, linearized in 6 3 , is 



A e a(A, a(/x 2 ), e) 



a(// 2 



a 2 {n 2y 



+ 



r /(A e ,a(/i 2 ),6) , (27) 



l- 7 (A e ,e)a(/i 2 ) 6 2 [1 - 7 (A% e)a(/i 2 )] 2 ' 

with 7(A e , e) = ^(A- - 1), /(A e , a, e) = 1 - A- - (l + ^) ln[l - 7 (A e , e)a]. We 
shall now justify the boundary condition in eq. (]22]) . To this end, note that for 
e < (d > 4), the dimensionally-continued running coupling vanishes at zero scale, 
a(0, a(A* 2 ), e) = 0, order by order in its perturbative expansion in the coupling a(fif) 
evaluated at any nonzero scale Dimensionally-continued radiative corrections to 
therefore vanish at Q 2 = 0. 

The functions Wf\ Wf } can be determined by choosing v — 1 in eq. (|25| ) and 
expanding the Sudakov equation to second order in the (i-dimensional running cou- 
pling 

d 



cr 



ln^(Ar,Q7//X// 2 ),e) = a(Q7//, a(^), e) W^{N, 1, e) 



<9Q 2 



+ a 2 (Q 2 / yU 2 ,a( / u 2 ),e)H/f (iV,l,e) 



(28) 
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in terms of a s (p 2 ), using (see eq. (p7|)) 

A e a(A, a(/i 2 ), e) = a{p 2 ) + -f{X e , e)a 2 {p 2 ) + £>(a 3 ) 



(29) 



The resulting one- and two-loop coefficients of the evolution kernel can then be 
derived from low order calculations of the correction factors via 



w£\N,l,e) = (Q 2 /^yQ 2 ^- 2 p { ;\N,Q 2 /^ 2 ,e) 



(30) 



and 



W$\N,l,e) = (Q 2 /p 2 ) 



2 / ,.2\2e 



\d\nQ 2 



pf(N,Q 2 /p 



2 I ..2 



6 ~ 



~ti\N,Q 2 /p 2 



d 



l((Q 2 /p 2 )% e) Q^pfiN, Q 2 /p 2 , e)\. ( 31 ) 



Of course, unlike for Drell-Yan |fLl]| , the second order corrections to Higgs produc- 
tion have not yet been calculated. The above results can easily be inverted to derive 
Wfjf (z, 1, e), if the functions p^' 2 ^ are known. The full functions (z, £ 2 //i 2 , e) 

may then be constructed by reexpanding the running coupling a(£, 2 /p 2 ,a(p 2 ),e) in 
a(p 2 ), using eq. (|2T|) with A = £ 2 / 'p 2 , 



W£\z,e/» 2 ,e) 



^jW^iz^e), 

,2\ 2e 



(32) 



wfW//* 2 ,*) = fe) e < ) (z,i, e )+ 7 ((e 2 //i 2 r,6)(^) e < ) (z,i,6) 



In this way the one-loop coefficients W| can be straightforwardly determined 



from eqs. ([T3,|TP|). They may be written as 



W?\z, l,e) = 8(l- z)jf\e) + z* Lf}*yL ) + h (1 \z, e) 



(33) 



where the coefficient functions are regular functions of their arguments 

at z — 1. In the present case they are given by 



(i) ( 

,(1) 



/i (1) (^,e) 



_ /ll 203 7T 2 

(i) 



/TO = f ( HMe)-2C A e 



C A (l + ** + (1 - z)^ 
o 



(34) 
(35) 
(36) 

(37) 
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Given the limitations of the factorization theorem near the edge of phase space |p2 
from which the evolution equation (|2lD is derived [Q, we might be tempted to im- 
mediately discard all terms of order 1/JV. However, we wish to be careful with these 
terms, and we will examine them and their relevance later in this section. 

As was shown in WM, the term 5(1 — z) and the plus-distributions in eq. fl33|) 



are separately renormalization group invariant. We are therefore free to choose differ- 
ent functions v(z) for these two terms in the general resummed expression eq. (f25|). 
The natural choice for the 5(1 — z) n term is v — 1. Changes in fj, generate terms 
62 ln(//'/ //) at higher orders. The A integral may then be carried out explicitly for 
this term in eq. (|25|). For the plus-distribution term however, the natural choice is 
v(z) = (1 —z) 2 , just as for Drell-Yan. Then, using £ = (1 — z)/i in eq. fl32|) , we absorb 
the factor (1 — z)~ 2e in (|33|) into the boundary of the A integral in eq. (p5| ), involving 
only the running coupling. The term hS 1 ' is of (9(l/iV 4 ) and will be neglected. 

We now turn our attention to the function g^ l > (z, e) in (0). We have some freedom 
in its treatment, as different choices will only differ by constants or in OQ11N/N). 
However, we will argue that among all such terms the ones generated by g^(z, e) are 
universal and can legitimately be included in the resummed expression. To exhibit 
the importance of the different treatments of the residue function g^\z, e) we choose 
three schemes which probe the full range of possibilities. After rescaling to incorporate 
the factor (1 — z)~ 2t and combining the plus distribution with the Mellin transform in 
eq. fl2"5|), we see that the relevant function to approximate is (z N ~ x — l)g^(z,e) [We 



are neglecting the coefficient z e , because it is contributing at 0(1/N).}. The three 
schemes are defined by 

scheme a : -^(z N ~ l - l)g {1) (z, t) -> (z N ~ l -1)2 
Ca 

scheme (3 : ^r(z N ' 1 - l)g {1) (z, t) -> (z N ~ l - 1) 2 - (1 - z)(2z 2 - Az - 2z 3 ) 

scheme 7 : -^r(z N ~ l - l)g (1) (z, e) -> {z N ~ x - 1) 2 - (1 - z)(2z 2 - Az - 2z 3 ) 

-Az N - l (l-z). (38) 

The minimal scheme a involves replacing g^(z, e) simply by g^(l,e), scheme (3 
includes all terms of 0(1) in the exponent, whereas scheme 7 includes in addition 
some 0(ln*A^/A^) terms in the exponent. Using the one-loop evolution kernel W9~ 
we can now construct the resummed expressions for the Higgs production correction 
factor in the three schemes. However, these expressions are still divergent for e — » 0. 
The divergences are cancelled by mass factorization and renormalization for which 
we choose the MS scheme |23[ throughout this paper. For this purpose we need the 
resummed MS gluon distribution El. To one loop accuracy it may be written as 
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14 



24 



MS 



A/, — ,a(/i ),e 



cxp 



/u 2 



/2 



/X 

M 2 
A* 5 



r$(tf,a.(/i' 2 )) 



/2 99 



M 2 



where 



= exp 
x exp 



2 Jo 



dz 



z 



N-l 



- 1 



1-z 

f2 I, ,2 



^ 2 /m 2 dA 



(39) 



o 



A 



a(A, a(/i 2 ), e) 



and 



20 



11 CU /-^dA a 
/ M 2 



, a(/z 2 ),e 



exp 



a^/m 2 dA 

T 



a(A, a(yU 2 ), e) 



(40) 



(41) 



where M is the mass factorization scale. As is well-known, the one-loop anomalous 
dimension Vyj(N,a a ) is derived from the residue of the collinear singularity in the 
gluon operator matrix element. The function zq is related to that component of the 
residue which is proportional to the one loop coefficient of the QCD /^-function. The 
strong coupling in the LO cross section (|7]) has been left unrenormalized so far. The 
renormalization of this bare coupling is now performed in the MS scheme, via the 
replacement 



R 2 

a s ,B = a s (R 2 ) (i? 2 ) Z Q I — ,a s (n 2 ),e 



(42) 



a 



.{R 2 ) (R 2 ) 



exp 



T 



{b 2 a(\, a^j 2 ) , e) + b 3 a 2 (A, a(/i 2 ), e) + . . .} 



where we explicitly show the renormalization scale R. Thus the renormalization 
factorizes into several pieces according to the perturbative expansion of the QCD j3 
function, 



a s ,B = a. 



,{R 2 ) (R 2 )' 



z 0£ S (/J? 



R 2 



where 



Zi ^,a(/i 2 ),e 



exp 



« 2 /m 2 d\ 



y a(X,a(fi 2 ),e) 



(43) 



(44) 



Note that the above form of the Z-factor for the strong coupling constant, which 
is completely factorized from the correction factor, is very similar in form to the 
5(1 — z) piece of the MS density, and we may thus combine the overall coupling 
constant renormalization with the mass factorization procedure. Restricting ourselves 
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to next-to-leading order (i.e. putting z\ = 1) and choosing R = M one finds that 
simultaneous mass factorization and renormalization of the bare coupling in eq. ([?]) 
leads to 



P(P (N,Q 2 /M 2 ,a(M 2 )) 



0«(iV,MV^,a(^), e ; 



2 ' 



(45) 



Note that the dependence on the dimensional regularization scale /i drops out on the 
l.h.s. Because ps is infrared safe, we can return to four dimensions, and we find 



Ph,H (N,^,a(M 2 ) 



Pin (iV,^,a(M 2 ) 



Pin Uj^,*(M 2 ) 



exp 



-C* 



1 Z N-1 



dz- 



1 - z 



l 2 " 



x exp<^ a(Q 2 )C/ 



{i _ z)2 ^a(X,a(M 2 ),e) 

Q 2 



tT/3 + 203/36 - 11/6 In I 



-11/12 a 2 (Q 2 )C A b 2 In 2 



(iV,^,a(M 2 ) 



(46) 
(47) 



x exp 



-2C A ( dz(2z -z 2 + z 3 ) f 1 
Jo 7(i 



d-) 2 |k A 



a(\,a(M 2 



Pin U,j^,a(M 2 ) 



(4* 



x exp 



7o 



JV-1 



d-) 2 ^ "A" 



a(X,a{M 2 



frfn ( AT Q*_ a ( M 2^ x exp 



p,;& 7 at, 



M 2 ' 



2C>(Q 2 ) 



(49) 



A few remarks are in order. First, the above expressions are formally not well defined 
because the integration paths in the exponent traverse a singularity, related to an 
infrared renormalon. Our purpose in this paper is however not the numerical eval- 
uation of the resummed cross section; we consider the resummed formulae to be a 
generating functional for an approximation to the QCD perturbation series, rather 
than its approximate sum. 

The second remark is that one may try to incorporate (a part of) the two loop 
evolution kernel. Some of the NNLO terms in the exact calculation can be inferred 
using renormalization group methods, such as used for Drell-Yan in Ref. ]2"5| . A 
straightforward calculation along these lines leads to 



wi 2 \z,l,e) = z(l-z)~^P^(z) + 



Ml 



(50) 



where P^ is the two-loop component of the Altarelli-Parisi gluon-to-gluon splitting 
function |26|. Moreover, the mass factorization and renormalization of the overall 



12 



coupling in the Born cross section can be carried out in a similar fashion as at one 
loop, leading to the modified MS distribution (for R = M) 



Li Li Lb Lb 

(51) 

where 0^ denotes the resummed MS distribution based on the two-loop anomalous 
dimension T^j p6 |. However, we stress that in this paper we restrict ourselves to 
examining the consequences of resummation with the one- loop kernel W^-Q 

There are several arguments supporting the inclusion of In* N/N terms in the 
resummation via the 7-scheme. In momentum space such terms correspond either to 
the usual plus-distributions, which are already included, or to terms ln l (l — z), see 
the Appendix. The latter are usually discarded in Sudakov resummation procedures. 
Their inclusion extends the resummation formalism to incorporate subleading, in 
the sense that they are down by a power of iV compared to the plus distributions, 
divergent contributions for z — > 1 in addition to the plus distributions. 

Our first observation is that those ln l (l — z) terms which are included by using 
scheme 7 originate from the Altarelli-Parisi splitting function, as we have demon- 
strated above. In distinction to the plus-distributions, which have the same origin, 
they are not enhanced by the infrared eikonal factor 1/(1 — z), hence they are es- 
sentially of a collinear nature. Their origin thus implies that they are universal and 
independent of the hard process. A second observation is that the integral in the 
exponent in eq. (f£8"D, which generates these logarithms, is partly composed of the 
full MS distribution (fjfi|), which is an exact result, without any soft gluon approx- 
imation. A third supporting argument can be obtained from existing NNLO (MS 



scheme) calculations for the qq channel in the Drell-Yan reaction [|11|] and the non 



singlet structure function F2 |27[ i n deep-inelastic scattering. We are able to derive 
the leading ln l (l — z) terms and those subleading terms which are related to the 
running coupling at NNLO using the resummed formulae for the hard parts for these 



reactions, given in [L^,p^| (e.g. for Drell-Yan see eq. (5.18) and the analogous DIS 
result obtained from renormalizing eq. (5.13) in ||14|| ). 

Numerical consequences of including the ln J (l — z) terms will be presented in the 
next section for both Drell-Yan and Higgs boson production at the LHC. 



3. Two- and Three Loop Results 

In this section we examine the quality and the scheme dependence of our resummed 
formulae. We do this by expanding our results for the correction factors of Higgs boson 



^Except for a &2 ^(a^'/m) term at NNLO. 
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production and the analogous expression for the qq channel in Drell-Yan (eq. (5.18) 
in [14]), in various schemes (a, (3, 7) to NLO and NNLO in order to compare with the 
exact results of Refs. || and ||11|| . All comparisons will be made for the LHC We will 



find that for both reactions scheme 7 reproduces the exact results remarkably well, 
much better than the more conventional schemes a and (3 at low Higgs/7* masses. 
We further derive the scale-dependent logarithms in the NNLO expansion from our 
resummed expressions. Here also the agreement is reasonable. 

From the NNLO expansion of our resummed formulae for Higgs production, we 
can obtain the first indication on the importance of the NNLO contributions. As 
mentioned earlier, these estimates are important because the NLO K-factor is sizable. 
A full numerical study of the resummed cross sections, which requires a regularization 
of the renormalon singularity |28] , will be published elsewhere . 



Analogous to eq. flT5| ) the renormalized correction factor p<p of eqs. fl4T)| - f4"5D may 
be expanded as 



Q 2 /M 2 ,a) = £ a n (M 2 )pf\z,Q 2 /M 2 ) 



(52) 



n=0 



We begin with the NLO expansion. Using the Mellin transform formulae of the 
Appendix and the methods of Pfifl , we find that the 0(a) terms, upon inversion to 
momentum space [keeping the renormalization/factorization scale M 2 different from 
Q 2 ] are given by 



Pl$}(z, jtfi) = C A {m(z) + 2V (z)L M + (y + ^ - ^L M )5{1 - z)}, (53) 



,7T 2 203 11 



P?Sfe ^p) = C A AV ^ + 2T >o(z)L M + - z)} 



Pu 



^( Z ,^ 2 )=P^(z,^ 2 )-C A {8£ 1 ( Z )} 



M 2 



where we use the notation 



M 2 ' 

Q 2 



P$(z, Jfi) = P ( hkz, Jp) + 2«(1 - z) 



(54) 
(55) 

[for all three schemes] , (56) 



VAz) 



In* (1 



£ i (z) = W(l-z) 



ln(Q 2 /M 2 



(57) 



With the results above we construct the correction factors according to eq. ( [52] ) for 
the three schemes a, (3 and 7 so that we can compare with exact results. Similarly 
one may obtain results to 0(a 2 ). We find 



pI [ $^q 2 /m 2 ) 



+ 



C A {8C A V 3 (z) + (-46a + 12C A L M )V 2 (z 
(~^~c A — SC A ( 2 



4b?L 



22 

—C A L M + AC A L 



V x {z) 



14 



203 11 
+ (16CUCs - 4CUC 2 £m + —C A L M - b 2 L 2 M - -C A L 2 M )V (z) 

( 203 r t +Kr 1 or r± 41209 r 
^ 18~ 6C ^C 2 - 126X4 + 2592 C A 

203 11 2233 

+ 8wC3-^m ~~ ^C 2 b 2 L M — 7^rb 2 L M — —C A ( 2 L M — ttttC/J-'m 

3d 3 216 

+ ^C A L 2 M -2C A ( 2 L 2 M + ^b 2 L 2 M )5(l-z)}, (58) 

p?g(z, Q 2 /M 2 ) = C A {8C A V 3 (z) + (-46a + 12C A L M )^ 2 (^) 
+ (-8«2-46 2 L M + 4C A L 2 vf )I? 1 (^ 

+ (16^3-^2^-62^)^0(2) 

2909 

+ (6CX2 2 - 12C A C 4 - ^-r&2 - 2C 2 6 2 Lm 

216 

+ 8C A ( 3 L M -2C A ( 2 L 2 M )8(l-z)}, (59) 
pTJ(z,Q 2 /M 2 ) = p^(z,Q 2 /M 2 ) + C A {-WC A £ 3 (z) 
+ (862 + 8C A - 2AC A L M )£ 2 (z) 

+ {l§C A C, 2 + 8b 2 L M + 8C A L M -8C A L 2 M )£ 1 {z)}, (60) 

pftfaCf/M 2 ) = pt { S{z,Q 2 /M 2 ) + C A {%C A V l {z)+AC A V G {z)L M 

1 1 239 1 

+ (±C A ( 2 --C A L M + —C A -2b 2 L M )8(l-z)}, (61) 

pf\z,Q 2 /M 2 ) = p^(z,Q 2 /M 2 ) + Ca{8C a Pi(s) + 4C A V (z)L M 

+ (AC A ( 2 + 2C A -2b 2 L M )5(l-z)}, (62) 

pl (2) (^,Q 2 /M 2 ) = pf 2 I }(z,Q 2 /M 2 ) + C A {8C A V 1 (z)+AC A V (z)L M 

WC A £ 1 (z) + (4C A ( 2 + 2C A -2b 2 L M )5(l-z)}. (63) 

In constructing the scale logarithms in the above expressions we have used the re- 
placement 

a(Q 2 ) = a(M 2 ) - a(M 2 ) 2 b 2 L M (64) 
in the last exponents of eqs. (|46|j4^|) . We note that this effectively amounts to including 

(2) 

one term from the two- loop Sudakov evolution kernel . This term can be derived 
from the one-loop evolution kernel using the renormalization group, see the discussion 
below eq. fl3"7|). In Fig. [| we present the correction factors for SM Higgs production 
at the LHC, which coincide with the correction factors of MSSM scalar Higgs boson 
production for small tg/3, and in Fig. [| for MSSM pseudoscalar Higgs production, 
where we have identified M 2 = Q 2 . For all results in this section we used the CTEQ4M 
parton densities Ol, a two-loop running coupling constant, with = 202 MeV for 
the NLO and NNLO quantities and CTEQ4L densities with a LO strong coupling 
constant (Ajfo = 181 MeV) for the LO quantities. 

In Figs. f|a and [5]a the "partonic" K-factors, obtained from folding the correction 
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factors with NLO parton densities and using a NLO strong coupling for all orders 
of the cross sections, are presented. For comparison we show in Figs. [|b and f)]b 
the corresponding NLO "hadronic" K-factors [which include the small contributions 
from normalized to the LO cross sections evolved with LO parton densities and a s . 
Whereas the former indicate the effect of the higher order corrections to the partonic 
cross section, the latter exhibit the convergence of the perturbative approach to the 
physical (hadronic) quantities. We observe from Figs. |] and [5] that at NLO scheme 
7 reproduces the exact NLO calculation almost exactly for the full range of the SM 
Higgs mass M H > 65 GeV and the MSSM Higgs masses > 45 GeV, whereas the 
schemes a and (3 agree with the exact result only for 1 TeV (the agreement of 

scheme a in the intermediate Higgs mass range is accidental). Moreover, note that the 
NNLO corrections to the partonic cross sections in scheme 7 are still significant. Full 
NNLO predictions for hadronic cross sections require NNLO parton densities, which 
are not yet available. At NLO, a significant reduction of the hadronic K-factors 
compared to the partonic K-factors can be read off from Figs. |]b and |5|b, indicating a 
more reliable perturbative QCD expansion contrary to what Figs. |]a and |5]a suggest. 
We point out that should the size of the NNLO corrections to the physical cross 
section warrant concern about the convergence of the perturbative approach, our 
resummation method can provide a tool to control such large corrections. 

Besides the three-loop anomalous dimension, the determination of NNLO densi- 
ties requires NNLO calculations for the physical quantities included in a global fit. 



However, there are presently only a few exact NNLO calculations available [11,27,31 
An approximate way to proceed might be provided by the use of NNLO expansions 
of resummed cross sections, which at NLO have to approximate the exact results re- 
liably. For future high-energy hadron colliders we expect this to require the inclusion 
of the novel subleading contributions that have been discussed in our analysis. Once 
approximate NNLO results have been obtained for several processes, e.g. heavy fla- 
vor production at the LHC, and the three-loop anomalous dimensions for the NNLO 
evolution of the parton densities have been calculated, a global fit of approximate 
NNLO parton densities can be performed. The same procedure could of course also 
be followed at the resummed level. 

To investigate the reliability of scheme 7 at NNLO we confront in Fig. |^ the 
approximate partonic K-factor, using NLO parton densities and strong coupling in 
all expressions, for the Drell-Yan qq production channel at the LHC as a function of 



the off-shell photon mass Q with the exact calculation of Ref. [11]. We do this at 



NLO and NNLO in schemes f3 and 7. The expressions used are given in eqs. ( [A. 77 ) 



and ( |A.78| ) of the Appendix. We observe that for scheme 7 the agreement is again 



excellent at NLO. As for the Higgs case, this is perhaps not so surprising, since the 
NLO answer is to a large extent included in the evolution kernel. However, note that 
the agreement is remarkably good even at NNLO. Clearly the ln ? (l — z) terms are 
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Figure 4: a) Exact and approximate two- and three-loop partonic K-f actors, convoluted 
with the NLO gluon-gluon luminosity dL 9 ^ LO j dr , in the heavy top-mass limit. The 
results for the three different schemes are presented as a function of the scalar Higgs 
mass M H , using NLO CTEQ4M parton densities ^ and a s /A^L = 202 MeV]. b) 
Hadronic NLO K-factor using LO CTEQ4L parton densities / |73|/ and a s [A^o = 181 
MeV] for the LO cross section and including the NLO contributions from kh- 
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dominant in this order as well. In scheme a we found at NLO a similar accidental 
agreement with the exact result, as in the Higgs case. In Fig. ^ we show the NNLO 
result for scheme a, which fails to approximate the exact NNLO curve, in contrast 
with scheme 7. Note that scheme f3 fails at medium and low Q 2 . The impressive 
agreement in NNLO, in combination with the general arguments given at the end of 
the previous section, gives us confidence to assert that the NNLO expansion for Higgs 
production in scheme 7 will be close to the exact value. 
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Figure 6: Exact and approximate one- and two-loop partonic K-factors, using NLO 
CTEQ4M parton densities and strong coupling = 202 MeV] in all orders 
of the cross section, for Drell- Yan, in three different schemes, as a function of the 
7* mass Q. The lower solid line is the exact NLO result, in the qq channel, and the 
upper solid line is the NNLO one. 



Next we investigate the consequences of the scale logarithms for Higgs production. 
It was found in Ref. f| that the scale dependence of the NLO cross section is still a 
monotonous function of the scales. In view of the outstanding agreement of the 
exact results with our approximate ones in scheme 7 for M 2 = Q 2 for the r = Q 2 /S 
dependence, we may use the same results to examine the scale dependence at NNLO. 
In fact, from arguments such as given in p5 |, one may deduce that eqs. (|58| - |63|) 
approximate the exact scale dependent terms very well, the only term lacking being 
proportional to the two-loop anomalous dimension eq. (|50l), which we have omitted. 
Again we use the Drell- Yan case to gauge the quality of scheme 7 in describing the 



19 



scale dependence. This is done in Fig. [7|for two values of the 7* mass Q. To obtain the 
best approximation, we have defined the curve labelled 72, the NNLO approximation 
in scheme 7, by adding the NNLO term of the expanded resummed exponential to 
the exact NLO result, rather than the 71 curve. We see that at NLO and NNLO 
scheme 7 again describes the scale dependence quite reasonably. Again, for both the 
exact and approximate expressions a consistent analysis requires the use of NNLO 
parton densities in determining the scale dependence, but these are not available 
yet. (It was shown in [] 1 that for the full NNLO Drell-Yan cross section, including 
other production channels, there is a strong indication, albeit based on NLO parton 
densities, that the scale dependence is significantly reduced compared to NLO). 

In Fig. H we present the scale dependence for SM Higgs production. The curve 
labelled 71 includes the sum of the approximate term p^ (1) of eq. (|D and the NLO 



contribution to kh of eq. ( |13|) for the gg initial state. The approximate NNLO 
result, labelled 72, has been obtained by adding the NNLO term p 1 ^ of eq. ( j60| ) 
and the corresponding contributions up to NNLO to kh of eq. ([HJ) to the exact 
NLO result [this significantly improves the approximation similar to the Drell-Yan 
process]. The curves labelled NNLO include the full NNLO scale dependence in the 
partonic cross section, which has been obtained from the exact NLO result by means 
of renormalization group methods, neglecting quarks at all stages. This curve has 
been normalized to the 72 curve at £ = 1. We observe that the NLO term of scheme 
7 deviates from the exact NLO slightly at large scales and significantly for small 
scales. This is caused by terms of 0(1/ N), which have been neglected. At NNLO 
there is a strong indication for a significant stabilization of the theoretical prediction 
for the total Higgs production cross section at the LHC. There are deviations between 
the NNLO and 72 curves at small and large scales, which are again due to terms of 
0(1/N) that have been neglected in scheme 7. We found similar results to hold for 
the MSSM pseudoscalar Higgs case. 

Finally, let us comment on the phenomenological implications of our results for 
the Higgs K-factor at NNLO. When the NLO corrections to the Higgs production 
cross section were calculated both for the infinite mass limit || and for the general 
massive case |J it was found that the ratio of the NLO cross section to the LO 
one could be larger than two if one used NLO parton densities and strong coupling 
for both cross sections at LO and NLO. In order to estimate the increase in size of 
the QCD corrected physical cross section the hadronic K-factor has to be defined 
by including the corresponding cross sections evaluated with parton densities and 
strong coupling at the same order, i.e. LO cross sections with LO quantities and 
NLO with NLO quantities. This hadronic NLO K-factor amounts to 1.5-2.0 in the 
phenomenologically relevant Higgs mass range. This indicates that the procedure to 
predict with increasing accuracy the physical cross section, by including higher order 
corrections consistently in all quantities entering the factorization theorem, seems to 
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converge better than one might conclude from Figs. f|a and ||a. A consistent NNLO 
hadronic K-factor requires NNLO strong coupling and parton densities, which are not 
yet available, but one might expect a further reduction from them. However, as long 
as we do not know the parton densities beyond NLO, we observe from Figs. |]a and |5|a 
that NNLO corrections to the partonic cross sections are sizable. Our demonstration 
that the bulk of the corrections originates from soft and collinear gluon radiation, 
and the fact that they can be resummed analytically, provides then a different way to 
organize the perturbative expansion in the phenomenologically relevant Higgs mass 
range: one may redefine the original QCD perturbation series by rewriting it as 
the product of our resummed expression, times a new series. Due to the excellent 
approximation of the original series by the expanded resummed series, the new series 
is expected to be very well behaved perturbatively. This is of course the standard 
method for making sense of perturbative QCD near the edges of phase space, where 
QCD corrections are large. However, the result can now be extended much further 
away from threshold due to the inclusion of the novel subleading contributions. Note 
that this procedure would also require resummed parton densities. We furthermore 
note that the large size of the corrections compared to Drell-Yan is partly due to the 
fact that for every power of a s a color factor Ca = 3 appears. We have seen that for 
Drell-Yan in the MS scheme, where the corresponding color factor is Cf = 4/3 but the 
analytical structure of the soft gluon corrections is quite similar, the corrections are 
considerably smaller. The same phenomenon was observed for the case of resummed 
heavy quark production in Ref. [[32 1. These color factors are correctly included in the 



resummed formulae. 

As remarked earlier, the evaluation of the resummed series has its own subtleties, 
related to the appearance of the infrared renormalon. For its treatment there have 
recently been a number of proposals ||28|| , which we will not discuss here. We anticipate 
that the resummed Higgs production cross section will be quite sensitive to the details 



of handling the renormalon, due the large color factor in the exponent 32 . 



4. Conclusions 

In this paper we have performed the all order resummation of soft gluon effects 
in Higgs production both for the Standard Model and its minimal supersymmetric 
extension, to next-to-leading logarithmic accuracy. We have extended the usual re- 
summation formalism to include logarithms which, although integrable, diverge in 
the partonic cross section near the edge of phase space. By expanding our resummed 
results to NLO and NNLO, and using the Drell-Yan process for comparison, we have 
shown that this extension expands the applicability of resummation efforts into the 
phenomenologically relevant Higgs boson production range at the LHC. An accurate 
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assessment of the expected Higgs production rate is of paramount importance for the 
LHC physics programme. However, a physical prediction of the NNLO cross section 
requires knowledge of NNLO parton densities, which is not yet available. Clearly, in 
this regard, it would be interesting to investigate the applicability of our extended 
formalism to many other QCD production processes with potentially large K-factors 
at NLO, e.g. heavy quark production (33[ both at the Tevatron and the LHC, or to 



revisit the Drell-Yan process for phenomenological studies along the lines of [|34]j. 
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Appendix A: Useful formulae 



In this appendix we collect some useful formulae used in section 2. We 
extending the Mellin transform table of Ref. p5[ up to 0(1/N). Define 



1 — X 



For the lowest four values of n this integral is, up to 0(1/N) 
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with = Ne lE and 7# denoting the Euler constant. Define also 



J n (N) = dx x 
Jo 



N-l in 



ln n (l - x) 



For the lowest four values of n this integral is, up to 0(1/A^) 



UN) = - 
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(A.65) 

(A.66) 
(A.67) 
(A.68) 

(A.69) 

(A.70) 

(A.71) 
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In JV 

JiW = -— (A.72) 

■MN) = ^ + § (A.73) 

MN) . -2* (A,4) 

Next we present the NNLO perturbative expansions of the resummed hard part u q q 
of the MS Drell-Yan cross section in two schemes, defined in analogy to eq. (|38|) . 
The relevant function to approximate here is (z N ~ x — l)g^ Y (z, e) with g^y(z, e) = 
Cp (1 + z 2 ). The schemes a, (3 and 7 are defined by the replacements 

scheme a : — — {z N ~ x 
C F 

scheme (3 : — — {z N ~ l 
C F 

scheme 7 : —-(2 _1 
G F 

-2^- x (l-z). (A.75) 

We find the results 

u^(z,Q 2 /M 2 ) = 5(l-^) + a(M 2 )C F {4P 1 (^)+2P (z)L M 

+ (2C 2 - - z)} + a(M 2 ) 2 C F {8C F P 3 (^) + (-46a 

+ 12CVL A/ )X? 2 (z) + (-8CK2 - 2CV - 46 2 L M + AC F L 2 M )V x {z) 

+ (16C F Cs - 4C f C 2 ^a/ - b 2 L 2 M - C F L M )V (z) 





- (^ iV - 1 -l)2 






-> (z JV " 1 -l)2 + (l- 






-> (z 1 "- 1 -!) 2 + (l- 


-*)(! + * 



1 1 



M 



2 

+ 8C F ( 3 L M -2C F ( 2 L 2 M )8(l-z)} (A.76) 
ufs^Cf/M 2 ) = S(l-z) + a(M 2 )C F {4V l (z)+2V (z)L I 



+ (2( 2 -4+h M )5(l- 

+ a{M 2 ) 2 C F {8C F V 3 {z) + (-46a + 12C i? L M )P 2 (z) 

+ (-8CVC2 - 16C F + 6C F L M - Ab 2 L M + 4C F L 2 w )D 1 (z) 

+ (I6CVC3 - 4C f C 2 ^m - SC F L M - 6 2 L 2 7 + 3C F L 2 , / )P (^) 

15 

+ (6CK2 2 - 12C F C4 - y& 2 - 8CK2 + 8C F - QC F L M 
+ 4b 2 L M — 2( 2 b 2 L M + 3C F ( 2 L M + 8C F ( 3 L M 

-b 2 L 2 M - 2C f C 2 ^m + \c F L 2 M )5(l - z)} (A.77) 



4 

u^(z,Q 2 /M 2 ) = u p , qg -(z,Q 2 /M 2 ) + 
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+ a(M 2 )C F { -45i(z)} 

+ a(M 2 ) 2 C F { - 8C F £ 3 (z) + (46 2 + 8C F - 12C F L M )£ 2 (z) 

+ (16C i , + 8C F C2 + 2C F L M + 46 2 L M -4C F L 2 ^)£: 1 (z)}. (A.78) 
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